Extra condition is necessary to have a unique cluster wave vectors set in the periodic 

cluster approximations. 
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We added an extra condition, original lattice symmetry of chosen cluster around cluster central 
site, to the cluster approximation methods with periodic boundary condition such as dynamical 
cluster approximation (DCA), effective medium approximation (EMSCA) and nonlocal coherent 
potential approximation (NLCPA). For each cluster size, this condition leads to a unique cluster 
wave vectors set in the first Brillouin zone (FBZ) where they preserve full symmetry of first Brillouin 
zone around K = 0. In this case whole cluster wave vectors are restricted to the FBZ and when 
number of sites in the cluster is equal to the whole lattice sites, these approximations recover original 
lattice symmetry in real and k-spaces. 



I. INTRODUCTION 

Different single site approximations such as coherent 
potential approximation (CPA)i, cluster approximations 
such as cluster CPA with open boundary condition, dy- 
namical cluster approximation (DCA)2i2ii and effective 
medium super-cell approximation (EMSCA)^&£ with pe- 
riodic boundary condition are used to approximate self 
energy of interacting disordered systems. In the single 
site approximations the k dependent of self energy is ne- 
glected while in the cluster approximations such as DCA 2 
for interacting and disordered systems, NLCPA^ for dis- 
ordered systems and EMSCAi in general for interacting 
disordered systems the cluster wave vector, K, depen- 
dent of self energy are considered. The DCA originally 
constructed in the k-space by dividing the first Brillouin 
zone (FBZ) in to N c cells. The wave vectors at center 
of these cells are called cluster wave vectors and denoted 
by {Ki, Ktv c }- They claimed these cluster wave vec- 
tors, K, are corresponds to a real cluster with periodic 
boundary condition 2 .. But for some clusters there are 
many possible sets of cluster wave vectors, {K}^ where 
some of cluster wave vectors are in the higher Brillouin 
zones. Recently to consider contributions of such differ- 
ent sets of cluster wave vectors, {K}, the origin of FBZ 
fixed at K = (0, 0, 0) and a phase <j) added to the cluster 
orthogonality condition, 
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They claimed, when number of lattice sites N c — » oo the 
phase, <f>, goes to zero and this relation converts to the 
following original lattice orthogonality condition, 
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Their uncorrected assumptions and results are due to two 
factors, first their chosen cluster sizes haven't original 
lattice symmetry around cluster central site because this 



could lead to existence of different set of cluster wave vec- 
tors where some of wave vectors are in the higher BZs. 
Second, the Born von Karman periodic boundary condi- 
tion, k.L = 27m 9 , imply that maximum number of clus- 
ter sites, N c , should be the whole lattice sites N, not 
N c — > oo. So in this limit their defined <f> is not going 
to zero, hence their orthogonality condition Eq[T] is not 
converting to the original lattice orthogonality condition 
Eq|5] Since real lattice sites and allowed wave vectors 
in the FBZ have a one to one correspondence, hence for 
each cluster with N c sites, it must exist a unique set 
of N c cluster wave vectors K in the FBZ. Therefore we 
must add an extra condition to the periodicity of clusters 
with N c sites to have whole cluster wave vectors in the 
FBZ and just a unique set of super-cell wave vectors K. 
Although periodicity condition for such real space clus- 
ters are considered 2 -^^^ but symmetry of original lattice 
around central lattice site of clusters are not considered. 
This is the weakness point of DCA and NLCPA where 
such symmetry are not considered. By considering both 
periodicity and lattice symmetry conditions for clusters, 
size of super-cells limit to those clusters that we can de- 
fine the Wigner-Seitz cell (WSC) for their central sites 
and if we translate this WSC by cluster sites vectors set, 
{r sc }, where defined with respect to the central cluster 
site, they cover whole cluster with out overlapping. For 
a 2 dimensional square lattice (2d) this is illustrated in 
Figure [TJ This important condition lead to interesting 
result where for each cluster there is a unique super- 
cell wave vectors, {K}, with full FBZ symmetry around 
K = (0,0,0). Here we reformulate the EMSCA in the 
real space to clarify this problem. 



II. MODEL AND FORMALISM 

We start our investigation by a general tight binding 
model for a disorder alloy system which is given by, 



H = 
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+ ( £ * ~m)cL c «t> ( 3 ) 

i<j 

where cj CT (ci CT ) is the creation (annihilation) operator of 
an electron with spin a on lattice site i and hi a — c\ a Ci a 
is the number operator, tl? are the random hopping 
integrals between i and j lattice sites with spin a respec- 
tively. \i is the chemical potential and £j is the random 
on-site energy, where takes with probability 1 — c for 
the host sites and S with probability c for impurity sites. 

The equation of motion for electrons corresponding to 
the above Hamiltonian, Eq[31 is given by, 

( {E - Si + fi)6u - t a ) G(l, j) = S l3 (4) 

i 

where G(i,j) is the random single particle Green func- 
tion. EqfJ] could be rewrite as 

G(E; i,j) = G a (E; + ]T G°(E; i, l) £l G(E; (5) 

i 

The Dyson equation corresponding to Eq[4] for the exact 
averaged Green function, G(l,j), isiS, 

G(E; i, j) = G°(E; i, j)+£ G°(E; i, l)E(E; I, l')G(E; 

w 

(6) 

where the real space self energy matrix, S, is defined by, 

(eG) = SG (7) 

By eliminate clean system Green function G° between 
Eqs[H]and[71 the exact random Green function G could 
be expand in terms of exact average Green function G 
and exact self energy £ as, 

G = G + G(£-£)G. (8) 

Note that although Eqs 03 [5] and [5] are exact but no exact 
solutions for them. So they should be solved approx- 
imately. In the next section by adding original lattice 
symmetry condition to the chosen super-cell, we reformu- 
late the effective medium super-cell approximation (EM- 
SCA). In this new formalism we show for any super-cell 
size with periodicity and original lattice symmetry con- 
ditions not only there is a unique set of super-cell wave 
vectors K but also all super-cell wave vectors occur in 
the FBZ with its symmetry. 

III. REFORMULATION OF THE EFFECTIVE 
MEDIUM SUPER-CELL APPROXIMATION 

In the EMSCA the real lattice is divided in to the simi- 
lar super-cells such that each super-cell have whole origi- 
nal lattice symmetry around it's central site. This means 
that for each super-cell, if any one sit on it's central lat- 
tice site, exactly see whole lattice symmetry!. Figure 
Ufa) shows this approximation for a 2 dimension square 



lattice which is divided in to super-cells with nine sites, 
N sc = 9, the central site of each super-cell denoted by 
green color. Figure [Hb) illustrate one of these super- 
cells. Each of these super-cells have original lattice sym- 
metry around its own central site. The set of sites inside 
each super-cell denoted by {/}. Note that it is possi- 
ble to divide this lattice in to super-cells with out lattice 
symmetry, for example in Figure [He) and (d) this system 
divided in to square super-cells with four sites, N sc = 4, 
which super-cell hasn't original lattice symmetry around 
non of its own sites. Therefore this cluster can not give 
use a correct results. Since for each lattice vector there 
is a unique wave vector in the FBZ in the k-space, so 
for the set of super-cell vectors set {r sc } it must exist 
a unique set of super-cell wave vectors in the k-space in 
the FBZ. To see lattice symmetry effect on the super-cell 
wave vectors, we reexamine the super-cell wave vectors 
corresponding to FigfT] (a) and (c). Fig f2] shows unique 
super-cell wave vector sets corresponding to nine sites ap- 
proximation, N c = 9, of FigQ](a). While for the four sites 
super-cell which is not preserved lattice symmetry, Figfl] 
(c) , there are four different super-cell wave vector sets as 
illustrated in FigO The super-cell wave vectors of each 
set are not preserve FBZ symmetry. Hence we found just 
super-cells with original lattice symmetry around their 
own central sites preserve a one to one corresponding be- 
tween real space super-cell vectors and super-cell wave 
vectors in the k-space. Therefore in the next part of this 
section we reformulate the EMSCA for super-cells with 
lattice symmetry. 

Since the super-cells in the real lattice are similar, in- 
tra super-cell multiple scattering kept and inter super-cell 
multiple scattering are neglected, hence the super-cell self 
energy, E SC (J, J), obey periodicity condition of the orig- 
inal lattice super-cells^, 

S sc (r /J + Lr c ") = S sc (r/j) (9) 

where L™ n is a vector connects central lattice sites of 
mth and nth super-cells, 

i 

where m, and ni are integer numbers and Nl c is num- 
ber of lattice sites along the ith primitive vector an. 
So N l ac B.i is length of super-cell along the and for a 
3-dimension system, number of sites in a super-cell is 
N sc = Nj c N ac N ac . mi — ni = k is an integer number. By 
applying Eqsl9land[T0lto the following exact self energy 
relations, 

q 

we have 

J2 S sc (q)e tq ' (rjJ+L -" ) = ]T £ sc (q)e™. (12) 
q q 



3 



(a) 



71 




(b) 




© © 


© © 
© @ 


© © 




@ 


© @ 
@ @ 


© © 
© © 


• @ 

© 


@ 

A/B 


© © 


© © 

© © 














©) ^/B) 


(d) 



FIG. 1: (Color on line) (a) and (b) Show dividing a two di- 
mension square real lattice in to similar super-cells of N c — 9 
with original lattice symmetry around their own central sites. 
The green color sites illustrate central sites of super-cells and 
hatched squares show the Wigner-Seitz unite cell (WSUC) 
of each super-cell of this lattice. The WSUC and super- 
cells have whole original lattice symmetry around their ori- 
gins. The super-cell vectors are r sc = mae x + nae y where 
m,n — —1,0,1 and a is lattice constant. By translate each 
WSUC with respect to super-cell vectors r sc = mae x +nae y , 
whole super-cell is cover with out overlapping, (c) and (d) 
show this 2d lattice divided in to four site clusters, N c ~ 4, 
where non of these clusters preserve original lattice symmetry 
around their own sites. 
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FIG. 2: (Color on line) Show FBZ of a 2d square lattice and 
its cells corresponding to super-cell with N c — 9. The unique 
set of super-cell wave vectors are K mn = m^-e x + n-£-e v 
where m and n take -1, and 1. All these wave vectors are 
in the FBZ and the super-cell wave vectors, K, have full FBZ 
symmetry around origin K = 
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FIG. 3: (Color on line) Show four 

different super-cell wave vector sets 

{(0,0),(0,7r),(*r,0),(7r,7r)} ) {(0, 0), (0, tt), (-tt, 0), (-tt, tt)}, 
{(0,0),(7r,0),(0,-7r),(7r,-7r)} and 
{(0, 0), (— tt, 0), (— tt, — tt), (0, — 7r)} corresponding to four 
sites cluster approximation of Fig[T] (c). The red square 
illustrate FBZ. Non of these cluster wave vector sets are in 
the FBZ and the super-cell wave vectors, K, haven't FBZ 
symmetry around origin K = 0. a — 1 is the lattice constant. 



where Eqfl2l imply that, 

e <qX"» = 1 ( 13 ) 

hence, 

q.L™ = 2W (14) 
where m' is an integer number. Eq |14l for the nearest 
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neighbour super-cells, n^mi — 0,1, defines a unique set 
of N sc super-cell wave vectors in the first Brillouin zone 
(FBZ) which have FBZ symmetry around K = 0. To 
avoid confusing them with the lattice wave vectors k, 
they are denoted by K = {K n }, 

K.L SC = 2vrm' (15) 

where L sc = J^. n^^aj. From EqOj)]for each super-cell 
size we can obtain super-cell wave vector set {K„}. It 
should be emphasize that the super-cell wave vectors set, 
{K„}, around center of FBZ have original FBZ symmetry 
due to original lattice symmetry of super-cell. To reex- 
amine this method we consider a limit case, when each 
super-cell extended to the whole lattice, that is N l sc = N l , 
hence N sc — N and L sc = L. So the translation vector 
between nearest neighbour super-cells, where each super- 
cell are included whole lattice sites, are 

L = J2^N l ai . (16) 

i 

EciITSl by considering EqUH] converts to the following pe- 
riodicity condition 

K (JV , o ^ w) .L = 2W. (17) 

The super-cell wave vectors set K(jv sc ^jv) defined by 
EqfTTI are exactly the lattice vectors k defined by the 
following periodicity Born von Karman condition, 

k.L = 2W (18) 

hence in this limit K^^n) = k and EMSCA become 
exact. 

In the next part of this section by applying EMSCA to 
the exact system we drive relation between real space and 
K-space super-cell self energies, average Green functions 
and the orthogonality relation. First we investigate self 
energies. In the EMSCA^ just self energy inside of each 
super-cell sites are non zero, 

V(i,j) = 5 iI 8jjV sc {I,J) (19) 

where I and J are restricted to the same super-cell. We 
apply EqQj)] to the following exact relation to find rela- 
tion between K-space and real space super-cell self ener- 
gies, £ SC (K„) and E SC (I,J), 

s ( k ) = ^E s ( j ^> irsj " k - ( 2 °) 

y 

Since number of super-cells in the whole lattice is j^- and 
in the super-cell approximation just self energies inside 
of each super-cell, S SC (J, J), are nonzero hence by using 
EqsilandCESl EqdD] reduces to, 

E SC (K„) = J- S SC (I, J)e*»- K ". (21) 



By converting Eqj2T] to the real space we have 

^£ sc (K„)e- l V./- K " = (22) 

K„ 

SC 13 K„ 

where by consider the following super-cell orthogonality 
condition 

±-E eiTljKn=5lJ (23) 

Eql23l reduces to, 

£.c(J, J) = T^- ^(K n )e- r »' K " (24) 

Egs l!?l"ll2l{l andl2"4"lare a set of equations that can convert 
real space E SC (K„) and K-space super-cell self energies 
E SC (J, J). Although relation between real space exact self 
energy,E(i, j; E), and super-cell self energy, E sc (7, J; E), 
is known using Eg 1 191 but in the K and k-spaces it should 
clarify. To find relation between exact self energy E(k) 
and super-cell self energy E SC (K„) we defining k = K„ + 
k , where k restricted to a sub-cell in the FBZ symmetry 
around super-cell vector K„. Hence in following exact 
self energy relation, 

k 

the summation over k could be separated in to two parts, 
= 5Z k ' , where k' are wave vectors with respect 
to K„i hence, 

S ( r ^) = ^E e "" ,K "E E ( K « +k V rij - k ' ( 26 ) 

K„ k ' 

this exact self energy, £(/, J), could be convert to the 
super-cell self energy, E SC (I, J), if we assume 

e ik '- rjJ = l (27) 

and 

^E(K n + k') = -^-E SC (K„). (28) 

k' 

Hence we have 

E sc (rjj) = e ir/J ' K "S sc (K„) (29) 

SC K„ 

Note that similar to Eg s 1771 and |2"B1 are introduce in the 
DCA3^. Our next task is to find such relations for 
Green functions. 
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To find relation between average super-cell Green func- 
tion G(K ra ) and exact average Green function G(k), 



k 



where G(k) is 



G(k) = (G n " 1 (k)-S(k))- 



(30) 



(31) 



and Gq 1 (k) is bare system Green function. Now we write 
each wave vector k in terms of its nearest super-cell wave 
vector K„ and inside cell wave vector k', k=K n +k' hence 
J2 k = X1k„ Sk' ■ By using these and EqJ57| the exact 
average Green function G(I, J) in terms of k-space Green 
function G(k) is given by, 

G V> J ) = ^££^ k ) eM " e<r "" k ' (32) 



K„ k' 



by using Eq[27|the exact average Green function, G(I, J), 
converts to super-cell Green function, G SC (J, J), 

G. D (i, J) = lY, eiriJ ' Kn E G ( K - + k ') ( 33 ) 



where relation between K-space super-cell Green func- 
tion, G(K„), and k-space exact average Green function, 
G(k), is given by, 



VG(K„ + k') 



k 



N 

~n7c 



G sc (K n ). 



(34) 



similar to this equation derived in the DCA 
formalisrr&2.ii. So relation between K-space and 
real space super-cell Green functions is given by, 



G SC (I,J) 



1 



N 8C ^ 



G SC (K„) (35) 



To complete our set of equations for calculation self en- 
ergy and Green function, we apply the EMSCA formal- 
ism to the self energy definition equation, EqjTl where we 
take average over all super-cells except one, called impu- 
rity super-cell, hence the random potentials in all super- 
cells except impurity super-cell are replace by super-cell 
self energies^, 



S 5C S sc S s 

Sy v 



V 



(36) 



/ 



Figure Q] illustrate the EMSCA formalism for a Id alloy 
system. Therefore Eqj7]in the EMSCA reduces to, 



(e«Gjr> 



sc— imp 



(37) 
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FIG. 4: (color on line) Show effective medium super-cell ap- 
proximation method for a Id alloy system, (a) Real Id alloy 
system divided in to similar three sites super-cells, N c — 3 
where each super-cell have full original lattice. The green 
color sites illustrate central sites of super-cells. By taking 
average over all super-cells except one where called impu- 
rity super-cell we have one impurity super-cell in an effective 
medium, where medium constructed from replacing super-cell 
impurities by super-cell self energies, (b) A super-cell with 3 
sites, N c — 3, in an effective medium super-cell approxima- 
tion. 



By inserting Eq[37| in to EqJH] this equation reduces to 
the following relation for the super-cell impurity Green 
function for any arbitrary lattice sites i and j, 

G™p (i,j;E) = G sc (i,j;E) + 
J2 LL , G sc (i, L; E)(e L S LL , - £ SC (Z, L'; E)) x 
G™ p (L',j;E). 



(38) 



When i and j restrict to impurity super-cell sites, EqU 
reduces to, 

G™*> (I,J;E)=G SC (I,J;E) + 

J2ll> G sc (I, L; E){e l 8 w - E SC (L, L'\ E)) x 



where 



Gl™ p (L',J;E) 



(G i ™P(I,J;E)) = G sc (I,J;E). 



(39) 



(40) 



Although Eq[39] can be separated in to two equations 
by definition of cavity Green functior£ but for disor- 
der systems it is not necessary to defined such cavity 
Green function, although we must defined it in gen- 
eral for interacting disorder systems^. Eqs[34j [35l [33 
I3TJ1 and [3D] construct a complete set of equations to 
be solved self consistently to obtain super-cell self en- 
ergy ,E SC (L, L'; E), and average super-cell Green func- 
tion G SC (I, J; E). The algorithm for calculation of av- 
erage Green function G SC (I, J; E) in the EMSCA is as 
follows 

1- A guess for K-space self energies E sc (K n ) usually 
zero. 

2- Calculate the super-cell average K-space Green 
functions, G(K n ), by inserting E,5 C (K n ) in G(K n ) = 

Ek^GoHKn + kO-EfKn))-!. 
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FIG. 5: (color on line) Show comparison of average density 
of states of a Id alloy system for CPA, iV c = 1, EMSCA 
with super-cell sizes, N c = 3 and N c = 7 where super-cell 
have original lattice symmetry. The strength length, S — 4t, 
impurity concentration is c = 0.25 and /i = 0. The difference 
between DOS are due to nonlocal corrections. 

3- Fourier transform of K-space G sc (K n ) to obtained 
real space Green function G SC (I, J; E). 

4- Calculate Gl™ p (I,J;E) by inserting E SC (X, L'\ E) 
and G SC (I, J; E) in to EqEl 

5- Calculate average Green function G SC (I,J;E) by 
using Eg 1401 by taking average over all possible impurity 
configurations. 

6- Calculate new self energies S sc (i, L'; E) by inserting 
obtained G 4 ™ p (7, J; E) and G SC (I, J; E) in EqEH 

Now we apply our method to a 1 dimension (Id) alloy 
system. Since just super-cells with odd number of lattice 
sites preserve the lattice symmetry with respect to their 



middle site, we compared calculated average density of 
states for different super-cell size. Figure[S] illustrate av- 
erage density of states for super-cells iV c = 1, 3, 7. By 
increasing super-cell size more peaks in the average den- 
sity of states. The corrections are due to nonlocal ef- 
fects where neglected in the single site approximations, 
JV e = l. ' 



IV. CONCLUSION 

We have added an extra condition to the cluster ap- 
proximations such as Dynamical cluster approximation 
(DCA), non local coherent potential approximation (NL- 
CPA) and effective medium approximation (EMSCA), 
where the chosen cluster must have original lattice sym- 
metry with respect to its central site. This condition 
leads to a one to one correspondence between super-cell 
sites and super cell wave vectors in the first Brillouin zone 
hence a unique set of super-cell vectors for each cluster 
size. One of this vector is K = (0, 0, 0) and other super- 
cell vectors around this vector have full FBZ symmetry. 
While without this condition, for some cluster sizes, some 
of the cluster wave vectors locating in the higher Brillouin 
zones. This lead to two effects, first asymmetry of clus- 
ter wave vectors with respect to FBZ center hence when 
number of sites in the cluster goes to whole lattice sites 
the cluster wave vectors are not cover whole FBZ. Sec- 
ond, existence of many set of cluster wave vectors. In our 
formalism when number of cluster sites goes to number 
of lattice sites the super-cell vectors recover whole FBZ 
and EMSCA become exact. 
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